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Abstract 

The phase velocity of light is co-parallel to the direction of energy flow in classical vacuum. 
However, in certain uncommon materials, these two vectors can be oppositely directed, in 
which case the phase velocity is termed 'negative'. This negative phase velocity (NPV) gives 
rise to many exotic phenomenons, such as negative refraction, inverse Doppler shift and inverse 
Cerenkov radiation, and has technological allure. According to quantum electrodynamics, the 
presence of a magnetostatic field makes vacuum an anisotropic medium for the passage of light. 
Under the influence of a sufficiently strong magnetostatic field, vacuum supports NPV. Such 
ultrastrong magnetic fields are believed to arise due to dynamo action in newborn neutron 
stars and in binary neutron star mergers, for examples. In view of the possible occurrence of 
negative refraction, the influence of ultrastrong magnetostatic fields must be carefully taken 
into account in astronomical observations relating to neutron stars and associated gamma-ray 
bursts. 
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In the usual descriptions of light propagation, as learned by generations of students from stan- 
dard textbooks on optics [8] and electromagnetics [9], the phase velocity casts a positive projection 
onto the direction of energy flow, as provided by the time-averaged Poynting vector. That is, the 
phase velocity is positive. However, in certain special circumstances, it is possible for the phase ve- 
locity to cast a negative projection onto the time-averaged Poynting vector; i.e., the phase velocity 
can be negative. 

The transition from positive phase velocity (PPV) to negative phase velocity (NPV) is of partic- 
ular interest because NPV underpins the much-heralded phenomenon of negative refraction [1, 2], 
as well as many other exotic phcnomenons such as inverse Dopplcr shift and inverse Cerenkov radi- 
ation [3] . The scientific and technological possibilities offered by negative refraction such as the 
construction of near - perfect lenses from planar slabs of NPV-supporting materials have been 
widely reported on [2, 10]. Much of the research effort has been directed towards the realization 
of metamaterials which can be used to achieve negative refraction. These metamaterials are artifi- 
cial composite materials, often having complex micromorphologies [1, 11, 12]. Metamaterials with 
simple micromorphologies can also support NPV, as is demonstrated by certain homogenized com- 
posite materials [13]. A significant milestone was reached recently by experimentalists through the 
fabrication of metamaterials which support negative refraction at the optical wavelengths [11, 12]. 

We report here on the manifestation of NPV in a quite different context, namely in vacuum 
under the influence of a magnetostatic field = jS^] S^. In classical vacuum, the phase velocity 
is positive and the passage of light is unaffected by B^, as reported by an inertial observer. How- 
ever, this is not the case for the quantum electrodynamical (QED) vacuum. The QED vacuum is a 
nonlinear medium which can be linearized for rapidly time-varying plane waves. Thereby, for prop- 
agation of light, QED vacuum is represented by the anisotropic dielectric-magnetic constitutive 
relations [4] 



with eo = 8.854 x 10~^^ F m~^ and /xq = 47r x 10~^ H m~^ being the permittivity and permeability 
of classical vacuum, respectively. The relative permittivity and relative permeability dyadics of 
QED vacuum have the uniaxial forms 



D = eoe' E. 



B = ^oA* • K, 



(1) 




(2) 
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respectively, where / is the 3x3 identity dyadic and the constant a = 6.623 x 10~^^ H~^kg~^ s^. 
The constitutive dyadics (2) were derived by Adler [4] from the Heisenberg-Euler effective La- 
grangrian of the electromagnetic field [14, 15]. 
Suppose that plane waves with field phasors 

Mr) = exp {ik • r) , H{r) = Hj^ exp (ik • r) (3) 

and wavevector k = \k\k propagate through the QED vacuum described by (1) and (2). For 
simplicity, let and k both be aligned with the Cartesian z axis. In this case there is only one 
wavenumber, namely |fc| = ui-y/eojlo (unlike the B_g ^ k case where there are two wavenumbers for 
each propagation direction, as is described in the Supplementary Information section). The 

time-averaged Poynting vector P and phase velocity are straightforwardly derived by combining 
(1). (2) and (3) with the Maxwell curl postulates [16]. Thus, the projection of the phase velocity 
onto the time-averaged Poynting vector emerges as 

v,'P=^^{l-^am^)\E^\\ (4) 

Clearly, the phase velocity is positive for \B_J\ < l/VSa but negative for \B_J > l/VSa. In fact, 
the same inequalities for PPV and NPV hold when the direction of plancwavc propagation is 
perpendicular to the direction of the magnetostatic field, as is described in the Supplementary 
Information section. 

Magnetic fields of magnitude greater than l/VSa = 1.374 x 10^^ Tesla are needed for QED 
vacuum to support NPV propagation. To find fields of this magnitude we turn to astrophysi- 
cal environments. Ultrastrong magnetic fields, developed by dynamo action, are associated with 
certain neutron stars. For example, fields of the order of 10^°-10^^ Tesla, have been estimated 
for newborn neutron stars, such as soft gamma repeaters [5, 6]. Considerably stronger fields, of 
the order of 10^^-10^* Tesla, are predicted to arise during the merger of a binary neutron star 
system [7]. Accordingly, NPV propagation may be expected to occur in these neutron-star en- 
vironments. Across boundaries between two regions, one of which experiences a magnetostatic 
field of magnitude l/\/8a < 1.374 x 10^^ Tesla and the other a magnetostatic field of magnitude 
l/\/8a > 1.374 X 10^^ Tesla, negative refraction of light will occur. 

Our results have far-reaching implications for observational astronomy and theoretical astro- 
physics. The possibility of negative refraction arising from ultrastrong magnetic fields should be 
taken into consideration in estimating astronomical positions, particularly if the light path from 
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the detector to the point of observation traverses regions in the vicinity of neutron stars, for ex- 
ample. Furthermore, negative refraction may well influence the propagation of gamma-ray bursts 
and other electromagnetic radiation emitted by neutron stars. 
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Supplementary Information 

In our Letter we demonstrate that electromagnetic plane waves can propagate in a quantum elec- 
trodynamical (QED) vacuum with negative phase velocity (NPV), provided that a sufficiently 
strong magnctostatic field is acting. Here we provide further details of the plancwave analysis. 
As in our Letter, suppose that plane waves with field phasors (3) propagate through the QED 
vacuum described by (1) and (2). The wavevector k = \k\k and corresponding phasor amplitude 
Ejj = {Eq, Eq,E^) (and similarly ILg) are straightforwardly deduced by combining the frequency- 
domain constitutive relations (1) with the source-free Maxwell curl postulates 



V X H{r) + iu>D{r) = 

V X E{r) - iLoB{r) = 



(5) 



For our purpose, it suffices to consider only two cases: (i) propagation parallel to and (ii) 
propagation perpendicular to B^. Without loss of generality, our coordinate system is oriented 
such that B, = z. 
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For the case k = z, the dispersion relation yields one wavenumber, namely = w^eo/xo- The 
electric field phasor lies in the Cartesian xy plane but is otherwise arbitrary. The scalar product 
of the phase velocity and time-averaged Poynting vector is given as 

2Zp • £ = ^ (1 - 8a|B,H {\ES\' + lE^of) ■ (6) 

The case where ^ = x is more complicated as there are two possible wavenumbers, namely 
\k\ = fci,2 where 



, /l-20a|Bj2 /l-8a|BJ2 

fci=a;^W _ - fc2 = a;Vi^W _ - (7) 



The wavenumber fci corresponds to a propagating mode for \Bg\ < l/\/20a or \Bg\ > l/VSa; 
otherwise it is an evanescent mode. Similarly, the wavenumber k2 corresponds to a propagating 
mode for \Bg\ < l/\/24a or \Bg\ > l/\/8a; otherwise it is an evanescent mode. The electric 
field phasor is directed along the Cartesian z axis for the wavenumber ki, whereas it is directed 
along the Cartesian y axis for the wavenumber k2. The scalar product of the phase velocity and 
time-averaged Poynting vector is given as 

Vp-P=\ ° • (8) 

[ ^(l-24a|S,|2)|i?„^|2 for \k\ = k^ 

Therefore, plane waves propagate with NPV in QED vacuum provided that \B_g\ > l/\/8a, 
regardless of whether the propagation direction is parallel or perpendicular to B^. 
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